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Abstract 

A non-equilibrium picture of thermodynamics is discussed at the 
apparent horizon of FRW universe in f(R,T) gravity, where R is the 
Ricci scalar and T is the trace of the energy-momentum tensor. We 
take two forms of the energy-momentum tensor of dark components 
and demonstrate that equilibrium description of thermodynamics is 
not achievable in both cases. We check the validity of the first and 
second law of thermodynamics in this scenario. It is shown that the 
Friedmann equations can be expressed in the form of first law of ther- 
modynamics ThdS' h + ThdjS' = —dE' + W'dV, where d 3 S' is the en- 
tropy production term. Finally, we conclude that the second law of 
thermodynamics holds both in phantom and non-phantom phases. 

Keywords: Modified Gravity; Dark Energy; Apparent Horizon; Thermody- 
namics. 
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1 Introduction 

Cosmic observations from anisotropy of the Cosmic Microwave Background 
(CMB) [TJ, supernova type la (SNela) [2], large scale structure [3], baryon 
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acoustic oscillations [3] and weak lensing [5] indicate that expansion of the 
universe is speeding up rather than decelerating. The present accelerated 
expansion is driven by gravitationally repulsive dominant energy component 
known as dark energy (DE). There are two representative directions to ad- 
dress the issue of cosmic acceleration. One is to introduce the "exotic energy 
component" in the context of General Relativity (GR). Several candidates 
have been proposed [H]- [XQ] in this perspective to explore the nature of DE. 
The other direction is to modify the Einstein Lagrangian i.e., modified grav- 
ity theory such as f(R) gravity [TT] . 

The discovery of black hole thermodynamics set up a significant connec- 
tion between gravity and thermodynamics [12]. The Hawking temperature 
T = H^-, where n sg is the surface gravity, and horizon entropy S = ^ obey 
the first law of thermodynamics [32]- [E]. Jacobson [15] showed that it is 
indeed possible to derive the Einstein field equations in Rindler spacetime by 
using the Clausius relation TdS = 5Q and proportionality of entropy to the 
horizon area. Here 5Q and T are the energy flux across the horizon and Un- 
ruh temperature respectively, viewed by an accelerated observer just inside 
the horizon. Frolov and Kofman [16] employed this approach to quasi-de 
Sitter geometry of inflationary universe with the relation —dE = TdS to 
calculate energy flux of slowly rolling background scalar field. 

Cai and Kim [17] derived the Freidmann equations of the FRW universe 
with any spatial curvature from the first law of thermodynamics for the 
entropy of the apparent horizon. Later, Akbar and Cai [IE] showed that the 
Friedmann equations in GR can be written in the form dE = TdS + WdV at 
the apparent horizon, where E = pV is the total energy inside the apparent 
horizon and W = \{p — p) is the work density. The connection between 
gravity and thermodynamics has been revealed in modified theories of gravity 
including Gauss-Bonnet gravity [18], Lovelock gravity [HJ [20], braneworld 
gravity [21], non-linear gravity [22]-[27] and scalar-tensor gravity [191 EEJ 
[29] . In f(R) gravity and scalar-tensor theory, non-equilibrium description of 
thermodynamics is required [22]- [29] so that the Clausius relation is modified 
in the form TdS = 5Q + dS. Here, dS is the additional entropy production 
term. 

Recently, Harko et al. [30] generalized f(R) gravity by introducing an ar- 
bitrary function of the Ricci scalar R and the trace of the energy-momentum 
tensor T. The dependence of T may be introduced by exotic imperfect fluids 
or quantum effects (conformal anomaly). As a result of coupling between 
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matter and geometry motion of test particles is nongeodesic and an extra 
acceleration is always present. In f(R,T) gravity, cosmic acceleration may 
result not only due to geometrical contribution to the total cosmic energy 
density but it also depends on matter contents. This theory can be applied 
to explore several issues of current interest and may lead to some major dif- 
ferences. Houndjo [21] developed the cosmological reconstruction of f(R,T) 
gravity for f(R, T) = fi(R) + /2CO and discussed transition of matter dom- 
inated phase to an acceleration phase. 

In a recent paper [26], Bamba and Geng investigated laws of thermody- 
namics in f(R) gravity. It is argued that equilibrium description exists in 
f{R) gravity. The equilibrium description of thermodynamics in modified 
gravitational theories is still under debate as various alternative treatments 
[32] have been proposed to reinterpret the non-equilibrium picture. These 
recent studies have motivated us to explore whether the equilibrium descrip- 
tion can be obtained in the framework of f(R, T) gravity. The study of con- 
nection between gravity and thermodynamics in f(R,T) may provide some 
specific results which would discriminate this theory from various theories of 
modified gravity. 

In this paper, we examine whether an equilibrium description of ther- 
modynamics is possible in such a modified theory of gravity. The horizon 
entropy is constructed from the first law of thermodynamics corresponding 
to the Friedmann equations. We explore the generalized second law of ther- 
modynamics (GSLT) and find out the necessary condition for its validity. 
The paper is organized as follows: In the next section, we review f(R,T) 
gravity and formulate the field equations of FRW universe. Section 3 in- 
vestigates the first and second laws of thermodynamics. In section 4, the 
Friedmann equations are reformulated by redefining the dark components to 
explore the possible change in thermodynamics. Finally, section 5 is devoted 
to the concluding remarks. 



2 f(R, T) Gravity 



The action of f(R, T) theory of gravity is given by [30] 




(1) 
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where C( maU er) determines matter contents of the universe. The energy- 
momentum tensor of matter is defined as 



^matter) _ 2 5(^C { 

matter) J /„<. 

= ~^=f for ' 



We assume that the matter Lagrangian density depends only on the metric 
tensor components g^ u so that 



rp(matter) _ r __ ^ ^ {matter) /„\ 

1 fiU y^ut- {matter) q ■ 1<JJ 



Variation of the action ([I]) with respect to the metric tensor yields the field 
equations of f(R,T) gravity as 

R pv f R (R, T) - ^ v f(R, T) + (g^D - V^ u )f R (R, T) 

= 8nGT^ att ^ - f T (R, T )T^ atter) ~ MR, T)Q^ (4) 

where V M is the covariant derivative associated with the Levi-Civita connec- 
tion of the metric and □ = V M V^. We denote f R (R,T) = df(R,T)/dR, 

f T (R,T) = df(R,T)/dT and = The choice of f(R,T)=f(R) 

results in the field equations of f(R) gravity 

The energy-momentum tensor of matter is defined as 

T {matter) = ^ + + (5) 

where is the four velocity of the fluid. If we take H( ma tter) = ~Pm, then 

= -2T^ ter) -Pm9,u. (6) 
Consequently, the field equations (jlj) lead to 

R,J R (R,T) - ^g, u f(R,T) + (g, u D - V,V u )f R (R,T) 

= 8nGT^ atter) + T^ atter) MR, T) + p m g,J T (R, T). (7) 

The field equations in f(R, T) gravity depend on a source term, representing 
the variation of the energy-momentum tensor of matter with respect to the 
metric. We consider only the non-relativistic matter (cold dark matter and 
baryons) with p m = 0, therefore the contribution of T comes only from 



M1/ becomes 
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ordinary matters. Thus, Eq.(J7]) can be written as an effective Einstein field 
equation of the form 



]7iG eff T^ atter) + r 



(d) 



(8) 



where 

is the effective gravitational matter dependent coupling in f(R,T) gravity 
and 



T ,(d) _ 1 
^ f R (R,T) 



- 9tMU (f(R, T) - Rf R (R, T)) + (V„V„ - g^D)f R (R, T) 



(9) 

is the energy-momentum tensor of dark components. Here, prime means 
non-equilibrium description of the field equations. 
The FRW universe is described by the metric 



ds 2 = h aP dx a dx p + f W 



(10) 



where f = a(t)r and x° = t, x 1 = r with the 2-dimensional metric h a p = 
diag(—l,a 2 /(l — kr 2 )). Here a(t) is the scale factor, k is the cosmic curvature 
and dQ 2 is the metric of 2-dimensional sphere with unit radius. In FRW 
background, the gravitational field equations are given by 



3 H 2 + 



87rG eff p m + 



R 



: {Rf R - f)-3H{Rf RR (11) 



HI ■ Ml 1 ■ 4 I = — 



+ Tf RT ] 
1 



-(Rf R -f) + 2H(Rf RR + T^) + Rf RR 



R 2 Irrr + 2RTf RRT + Tf RT + T 2 f RTT 



These can be rewritten as 



3 H 2 + 



-2 H 



brG eff (p m + f/ d ) 
lirG eff (p m + p' d + p' 



(12) 

(13) 
(14) 



where p' d and p' d are the energy density and pressure of dark components 



1 



-(Rf R -f)-3H(Rf RR + ff RT ) 



(15) 



Pd 



\{Rf R ~f) + 2H(Rf RR + ff RT ) + Rf RR + R 2 f RRR 



+ 2RTf RR x + Tf R x + T 2 f R TT 



(16) 



Here T = 1 + ^^p- ■ The equation of state (EoS) parameter of dark 
fluid u' d is obtained as {p' d = 0J d p' d ) 
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RfiiR + R 2 Irrr + ^RTf RR T + Tf RT + T 2 f RTT — H(Rf RR + Tf RT ) 



(17) 



\{Rf R - f)-?>H{Rf RR + Tf RT ) 

The semi-conservation equation of ordinary matter is given by 

p + 3Hp = q. (18) 



The energy-momentum tensor of dark components may satisfy the similar 
conservation laws 



p d 4- 3H(p d + p d ) = q d , 
pt t + 3H(p tot + p tot ) = q tot , 



(19) 
(20) 



where p to t = Pm + Pd, Ptot = Pd and q to t = q + qd is the total energy exchange 
term and q^ is the energy exchange term of dark components. Substituting 
Eqs. CC?]) and (TH1) in the above equation, we obtain 



qtot 



8nG 



:h 2 + 



T 



(21) 



The relation of energy exchange term in f(R) gravity can be recovered if 
J= = 1. In GR, q^t = for the choice f(R, T) = R. 



3 Laws of Thermodynamics 

In this section, we examine the validity of the first and second law of ther- 
modynamics in f(R,T) gravity for FRW universe. 
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3.1 First Law of Thermodynamics 

Here we investigate the validity of the first law of thermodynamics in f(R, T) 
gravity at the apparent horizon of FRW universe. The dynamical apparent 
horizon is determined by the relation h al3 d a fdpf = which leads to the radius 
of apparent horizon for FRW universe 

ta = * (22) 

'H 2 + -4 



The associated temperature of the apparent horizon is defined through the 
surface gravity n sg as 

T h = (23) 



2tt 



where n sg is given by [T7] 



2\/-h r A 2Hr y 



(2H 2 + H+—). (24) 



In GR, the horizon entropy is given by the Bekenstein-Hawking relation 
Sh = A/AG, where A = 4irr\ is the area of the apparent horizon [T2]-|14j. 
In the context of modified gravitational theories, Wald |34j proposed that 
entropy of black hole solutions with bifurcate Killing horizons is a Noether 
charge entropy. It depends on the variation of Lagrangian density of modified 
gravitational theories with respect to Riemann tensor. Wald entropy is equal 
to quarter of horizon area in units of effective gravitational coupling i.e, 
S' h = A/4G e ff (35]. In f(R,T) gravity, the Wald entropy is expressed as 

S' h = (25) 
Taking the time derivative of Eq. (l22|) and using (fl4|) . it follows that 

f R df A = ^Gf- 3 A (p' tot + p' tot )HTdt, (26) 

where p' tot = p' m + p' d and p' tot = p' d . dfA is the infinitesimal change in radius 
of the apparent horizon during a time interval dt. Using Eqs. ff25|) and (126]) . 
we obtain 

VdS' h = 47rfi( P ; oi + p' tot )Hdt + -^-df R + T -4^d (±) . (27) 



2nf A 2GF 2G V T 
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If we multiply both sides of this equation with a factor (1 — r ^j^Hr a\ it 
follows that 

T h dS' h = A^\{p' tot + p' tot )Hdt - 2^\{p' tot + p' m )dr A + nr *™* 
nr 2 A T h f R J ( 1 



Now, we define energy of the universe within the apparent horizon. The 
Misner-Sharp energy [36] is defined as E = |g which can be written in 
f(R, T) gravity as [22 I2S] 

E' = (29) 

ZLr eff 

In terms of volume V = ^r\, we obtain 



3V /„„ k 

~~2 I ~ " ytot 



E '=^krA H2+ ^ ]=Vp '" (30) 



which represents the total energy inside the sphere of radius ta- It is obvious 
that E' > 0, if G e ff = j^- > so that the effective gravitational coupling 
constant in f(R,T) gravity should be positive. It follows from Eqs. ffTB"]) and 
([30]) that 

.=*rj , _/ sua* , a^2J ->= • f Adf R , fU/ flj /' 1 



dtf = -47rrl(p^ + ^JM + + ^ + ^- j . (31) 

Using Eq.flSD in (J2E]), it follows that 



which includes the work density W = h(Ptot~Ptot) W?]- This can be rewritten 



2 

as 



(32) 
tten 

T h dS' h + T h d 3 S' h = -dE' + W'dV, (33) 



rA M , t ^r r ^ffn\ _ HE' + S' h T h )Jf R 



where 

d A = -^rd + *™U«I ^ j = -^f^ ^ j • 04) 

When we compare the cosmo logical setup of f(R, T) gravity with GR, Gauss- 
Bonnet gravity and Lovelock gravity [18j-[20j, we obtain an auxiliary term 



8 



in the first law of thermodynamics. This additional term d 3 S' h may be in- 
terpreted as entropy production term developed due to the non-equilibrium 
framework in f(R, T) gravity. This result corresponds to the first law of ther- 
modynamics in non-equilibrium description of f(R) gravity [26] for f(R, T) = 
f(R). If we assume f(R,T) = R, then the traditional first law of thermody- 
namics in GR can be achieved. 

3.2 Generalized Second Law of Thermodynamics 

Recently, the GSLT has been studied in the context of modified gravitational 
theories [25]- [28]. It may be interesting to investigate its validity in f(R,T) 
gravity. For this purpose, we have to show that [28] 

S' h + d 3 S' h + S' ot > 0, (35) 

where S' h is the horizon entropy in f(R,T) gravity, djS' h = d t (djS' h ) and S' tot 
is the entropy due to all the matter and energy sources inside the horizon. 
The Gibb's equation including all matter and energy fluid is given by [38] 

T tot dS' tot = d(p' tot V) + p' tot dV, (36) 

where T to t is the temperature of total energy inside the horizon. We assume 
that T to t is proportional to the temperature of apparent horizon [25], [28], i.e., 
T tot = bT h , where < b < 1 to ensure that temperature being positive and 
smaller than the horizon temperature. 

Substituting Eqs. (l33|) and ( 136]) in Eq.( ]35]) . we obtain 

247T- 

S' h + d 3 S' h + S' tot = — > 0, (37) 

where 

E = (l-b)p' tot V + (l- b -)(p' tot +p' tot )V 

is the universal condition to protect the GSLT in modified gravitational the- 
ories [2H]. Using Eqs. (!T3~]) and ( TH1) . condition ([35]) is reduced to 



where 



X = 2(1 -b)H(H- ^)(H 2 + ^)f R + (2-b)H(H-^) 2 f R 
+ (l-6)(F 2 + i) 2 ^(&. 

Lb J~ 

Thus the condition to satisfy the GSLT is equivalent to X > 0. In flat FRW 
universe, the GSLT is valid with the constraints d t (y-) > 0, H > and 
H > 0. Also, J 7 and /r are positive in order to keep i£ > 0. If b = 1, i.e., 
temperature between outside and inside the horizon remains the same then 
the GSLT is valid only if 

For Eq.fllOD, the effective EoS is defined as u eff = -1-2(H - ^)/3(H 2 + £). 
Here H < -\ corresponds to quintessence region while H > represents the 
phantom phase of the universe. It follows that GSLT in f(R,T) gravity is 
satisfied in both phantom and non-phantom phases. This result is compatible 
with [39] according to which entropy may be positive even at the phantom 
era. Bamba and Geng [261 EZ] also shown that second law of thermodynamics 
can be satisfied in f(R) and f(T) theories of gravity. 



4 Redefining the Dark Components 

In previous section, we have seen that an additional entropy term d^S^ is 
produced in laws of thermodynamics. This can be considered as the result 
of non-equilibrium description of the field equations. If we redefine the dark 
components so that the extra entropy production term is vanished, then such 
formulation is referred as an equilibrium description. It has been seen so far 
that the equilibrium description does exist in modified theories of gravity 
[26l [271 [29] and extra entropy production term can be removed. 

Here, we discuss whether the equilibrium description of f(R,T) gravity 
can be anticipated. In fact, we may reduce the entropy production term 
through this description but it cannot be wiped out entirely. We redefine 
the energy density and pressure of dark components. The (00) and (11) 
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components of the field equations can be rewritten as 

k 



3 H z + 



2 H 



a" 
k 



8nG eff (p m + p d ), 
8nG eff (p rn + pd + Pd), 



(40) 
(41) 



where G 



eff 



G 



MR,T) 

8vr 



is the effective gravitational coupling, p d and pd 



are the energy density and pressure of dark components given by 



Pd 



+ 



SnGT 
k 

~2 



-(Rf R - f) - 3H(Rf RR + ff RT ) + 3(1 - f R )(H 2 



a 



(42) 



Pd 



8nGF 



" /) + 2H(Rf RR + ff RT ) + Rf RR + R 2 f RRR 



+ 2Rff RRT + ff RT + f 2 f RTT - (1 - f R ) (2H + 3H 2 



k 



(43) 



The EoS parameter uj d in this description turns out to be 

^d = — 1 + {Rf RR + R 2 f RRR + 2RT f RR x + Tf R x + T 2 fuTT — H(Rf RR 
+ T/rt) - 2(1 - f R )(H - ^)}/{l(Rf R -f)- 3H(Rf RR + ff RT ) 



+ 3(l-f R )(H 2 + - 2 )}, 

In this case the expression of total energy exchange is given by 

3 / TT o k 



Qtot 



8nG 



H 2 + -)d t 



(44) 



(45) 



Since dt(fr(R, T)) ^ in f(R, T) gravity, so that q to t does not vanish. So, we 
may not establish the equilibrium picture of thermodynamics in this modified 
gravity. Hence, again we need to consider the non-equilibrium treatment of 
thermodynamics. This result differ from other modified gravitational theories 
due to the matter dependence of the Lagrangian density. In f(R) gravity 
the redefinition of dark components result in local conservation of energy 
momentum tensor of dark components It is clear from Eqs.(ll7j) and 
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(|44p that the EoS parameter of dark components is not unique in both cases. 
Thus, one should consider both formulations of the field equations in cosmic 
discussions. 

Now we check the validity of the first and second laws of thermodynamics 
in this scenario. 



4.1 First Law of Thermodynamics 

In this representation of the field equations, the time derivative of radius r A 
at the apparent horizon is given by 

dr A = A-Kf\GF{ptot + PtotjHdt. (46) 

Since in f(R,T) gravity, the equilibrium description is not feasible as it can 
be seen in modified gravitational theories such that f{R), f(T) and scalar 
tensor gravity etc. Thus, we use the Wald entropy relation Sh = A/(4G e ff) 
rather than introducing Bekenstein-Hawking entropy. Using Eq. fj46|) . the 
horizon entropy becomes 

^3/ . \tt,. . r A 



—dS h = 4nf A (ptot + Ptot)Hdt + —d I — I . (47) 



2nr A 2G \T 

The associated temperature of the apparent horizon is 

r » = ^ 1 -^- (48) 

Equations fT4"7]) and fT4"gj) imply that 

7rf 2 Th ( 1 \ 

T h dS h = Aixf\( Ptot + p tot )Hdt - 2irr 2 A ( Ptot + p to t)dr A + [jj- ( 49 ) 

Introducing the Misner- Sharp energy 

E = ih = vptou (50) 

we obtain 

-^3/ , „ \ 7LT-7J. I „ ~2 _ ,~ , f A J 1 



dE = -A7Tf A (p tot + Ptot )Hdt + A7rr 2 A p tot dr A + ^d ^ - ) . ( r ) I ) 
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The total work density is defined as [37] 

W = -lr^h a0 = l -{ Ptot - Vtot \ (52) 

By combining Eqs. (fl9l) . (I5TI) and (|52|) . we obtain the following expression of 
first law of thermodynamics 

T h dS h + T h d 3 S h = -dE + WdV, (53) 

where 

^ = + 2 ^ d (j) = - T {h s ") d ( y) 

n{4H 2 + H + 3k/a 2 ) f 1 

~G(H 2 + k/a 2 )(2H 2 + H + k/a 2 ) 

is the additional term of entropy produced due to the matter contents of 
the universe. It involves derivative of f(R,T) with respect to the trace of 
the energy-momentum tensor. Notice that the first law of thermodynamics 
T^dSh = —dE + WdV holds at the apparent horizon of FRW universe in 
equilibrium description of modified theories of gravity [26j [2TJ, [29] . However, 
in f(R, T) gravity, this law does not hold due to the presence of an additional 
term d^Sh- This term vanishes if we take f(R,T) = f(R) which leads to the 
equilibrium description of thermodynamics in f(R) gravity. 




4.2 Generalized Second Law of Thermodynamics 

To establish the GSLT in this formulation of f(R, T) gravity, we consider the 
Gibbs equation in terms of all matter field and energy contents 

TtotdS tot = d{p tot V) + PtotdV, (55) 

where T to t denotes the temperature of total energy inside the horizon and 
Stot is the entropy of all the matter and energy sources inside the horizon. 
In this case, the GSLT can be expressed as 

S h + d 3 S h + S tot > (56) 

which implies that 

> 0, (57) 



bRGF(H 2 + £ 
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where 



y 



+ (i 



2(1 



- b)H(H - + A) + (2 - 6)tf(# - A) 



2 



Thus the GSLT is satisfied only if y > 0. In case of flat FRW universe, the 
GSLT is met with the conditions d t {y) > 0, H > and H > 0. In thermal 
equilibrium 6 = 1, the above condition is reduced to the following form 



for V = ~-kt\ and R = 6 (if + 2H 2 + k/a 2 ). <B > clearly holds when 
the Hubble parameter and scalar curvature have same signatures. It can 
be seen that main difference of results of f(R,T) gravity with f(R) gravity 
is the term T = 1 + ^^p- ■ We remark that in both definitions of dark 
components, the GSLT is valid both in phantom and non-phantom phases of 
the universe. 

5 Concluding Remarks 

The fact, f(R,T) gravity is the generalization of f(R) gravity is based on 
coupling between matter and geometry [301 • This theory can be applied to 
explore several issues of current interest in cosmology and astrophysics. We 
have discussed the laws of thermodynamic at the apparent horizon of FRW 
spacetime in this modified gravity. Akbar and Cai |23| have shown that 
the Friedmann equations for f(R) gravity can be written into a form of the 
first law of thermodynamics, dE = TdS + WdV + TdS, where dS is the 
additional entropy term due to non-equilibrium thermodynamics. Bamba 
and Geng [26l [27] established the first and second laws of thermodynamics 
at the apparent horizon of FRW universe with both non-equilibrium and 
equilibrium descriptions. 

We have found that the picture of equilibrium thermodynamics is not 
feasible in f(R, T) gravity even if we specify the energy density and pressure 
of dark components (see Section 4). Thus the non-equilibrium treatment 




(58) 
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is used to study the laws of thermodynamics in both forms of the energy- 
momentum tensor of dark components. In f(R,T) gravity, q to t does not 
vanish so there exists some energy exchange with the horizon. The non- 
equilibrium description can be interpreted as due to some energy flow between 
inside and outside the apparent horizon. The first law of thermodynamics is 
obtained at the apparent horizon in FRW background for f(R,T) gravity. 

We observe that the additional entropy term is produced as compared 
to GR, Gauss-Bonnet gravity (18] , Lovelock gravity [19j [20] and braneworld 
gravity |21j. The equilibrium and non-equilibrium description of thermody- 
namics in f(R) gravity can be obtained if the term Jt(R,T) vanishes i.e., 
Lagrangian density depends only on geometric part. We have established the 
GSLT with the assumption that the total temperature inside the horizon T tot 
is proportional to the temperature of the apparent horizon Th and evaluated 
its validity conditions. The GSLT in f(R) gravity follows from condition 
f )38|) if T = 1. It is concluded that in thermal equilibrium, GSLT is satisfied 
in both phantom and non-phantom phases. 
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